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Abstract 

We suggest an approach to the problem of finding integral equations 
for the excited states of an integrable model, starting from the Thermo- 
dynamic Bethe Ansatz equations for its ground state. The idea relies 
on analytic continuation through complex values of the coupling con- 
stant, and an analysis of the monodromies that the equations and their 
solutions undergo. For the scaling Lee- Yang model, we find equations 
in this way for the one- and two- particle states in the spin-zero sector, 
and suggest various generalisations. Numerical results show excellent 
agreement with the truncated conformal space approach, and we also 
treat some of the ultraviolet and infrared asymptotics analytically. 
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1 Introduction 



The Thermodynamic Bethe Ansatz (TBA) |T| has proved to be a very useful 
tool in the study of integrable two-dimensional field theories. The finite- 
volume ground state energy can now be found for many models, expressed 
in each case in terms of the solution of a set of nonlinear integral equations. 
These turn out to be vulnerable to both numerical and analytical attack, and 
yield a large amount of nontrivial information. Given these successes, it is 
natural to hope that similar equations might describe the remaining energy 
levels. However, the derivation of the TBA equations relies on arguments 
which perforce single out the ground state, and, save for a few states which 
become degenerate with the ground state in large volumes ||, the 
desired generalisation has proved elusive. 

In this paper we propose to sidestep the problem by returning to the old 
idea that one can move between energy levels by analytic continuation in a 
suitable parameter. In 0+1 dimensions, the quantum-mechanical problem, 
this was seen most spectacularly in the analysis of the quantum anharmonic 
oscillator performed by Bender and Wu |J. Somewhat simpler in analytic 
structure, but more relevant to the following, is the finite volume spectrum 
of the Ising field theory in 1+1 dimensions. If the 'spatial' dimension is 
rolled up into a circle of circumference R, then the ground state energy (the 
lowest eigenvalue of the infinitesimal transfer matrix in the 'time' direction) 
can be written as Q 

4 M (M, R) = EZ k (M, R) - ^(MR) , (1.1) 

where M the inverse of the bulk correlation length, E™ k (Al,R) contains 
the (as it happens logarithmically divergent) bulk contribution, and 

tm / \ 1 11 

4 W = log- + - + lnff-7E 




r 2 



r>U2k-l)V-(2k- l) ,- w —^- (1.2) 



with r = MR and 7^ = 0.57721566... Apart from the logarithmic singularity 
at R = 0, this exhibits a series of square root branch cuts, evenly spaced 
along the imaginary axis. Suppose that R is continued into the complex 
plane along a path enclosing the singularities at k = k\ , &2 • • • k n . Then on 
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its return to the real axis, Eq (M, R) has been replaced by 



2 n / 

E^ k2 _.. kn (M,R) = E™(M,R) + -^2^ + (2k-l)^ , (1.3) 

i=l 

an excited state with ultraviolet scaling dimension J2i(^i~ !)• This covers 
all of the symmetrical descendants of the primary fields I and e in the 
spin-zero sector. To find descendants of the spin field a, the same process 
can be repeated, but this time starting from the lowest excited state in the 
low-temperature regime, which degenerates with Eq in infinite volume and 
is accessible as the ground state with twisted boundary conditions (see for 
example Q). 

In more general situations explicit expressions such as flyD, O 

are 

not available. Even for integrable perturbations of conformal field theories, 
the best one can do is to express the ground state energy in terms of the 
solutions e a {9) to a set of TBA equations. The simplest case is the scaling 
Lee- Yang model, or SLYM. This field theory, a perturbation of the non- 
unitary minimal model A4(2/5) by its unique relevant operator <p, has the 
action 

-4slym = Am(2/5) +l ^J V{x)d 2 x . (1.4) 

With M (A) = (2.642944 . . .)A 5 / 12 § and r = M(X)R, the single TBA equa- 
tion reads |J 

e{9) = r cosh 9 -$*L(9) , (1.5) 

where 

L[0) = log (l + e~^) , f*g(6) = — J_J0'f{0 ~ 0')g(9) , (1.6) 

and 

a/o\ 9 i Qfa\ Gfa\ sinh(6>) + z sin(7r/3) 

= - l aa lo sS{d) , S(6) = , — . (1.7) 

09 smh(6>) - £Sin(7T/3) 

The function S{9) is the S'-matrix of the single neutral particle in the 
model Q, and has the '(/> 3 ' bootstrap property, that S(9—m /3)S(9+m /3) = 
S(9). In terms of these quantities, the ground state energy is 

E (X, R) = £ bu ik(A, R) - 7^co(r) (1.8) 

bit 



with 

£buik(A, R) = R , c (r) = ^J C °d9r cosh 9L{9) . (1.9) 



2 



At first sight (p~T5|) — ( JT79| ) are very different from Ql.l| ), (p~^), and it is not clear 
that analytic continuation will be practicable. To get a clue as to how to 
proceed, return to (|b^) and consider its alternative integral representation: 

c l M (r) = y°°^r cosh log (l+e- rcoshe ) . (1.10) 

As r moves into the complex plane, a singularity in c IM might be expected 
whenever i +e ~r^he = Q for 

some real 9q. However, deforming the contour 
of integration away from the real axis near 6q shows that such impressions 
are generally deceptive. This manoeuvre would only fail if two singularities 
were to approach the real axis from opposite sides, trapping the contour. 
This gives rise to a so-called 'pinch singularity', usually a branch point. 
If r is continued along some path encircling the critical value, the two 9- 
singularities (in this case, singularities in log(l+e _rcosh6 ')) execute a little 
dance in the complex plane (here, they just swap over), after which the 
contour has become tangled up. Undoing the tangle gives rise to the discon- 
tinuity across the cut. All of this is very clearly explained in chapter 2 of the 
book by Eden et al ||, and here we simply remark that it is an instructive 
exercise to recover the results already quoted by the use of such methods. 

For more general TBA systems, life is complicated by the replacement 
of rcosh# by e(9), a function which may itself be subject to nontrivial 
monodromies. Furthermore, e{6) is not known explicitly even for the ground 
state. Fortunately, only qualitative, 'topological', information about the 
movement of singularities is needed in order to deduce the modified TBA 
equations, and for this numerical work will suffice. 



2 One-particle states 

Both the truncated conformal space approach (TCSA) |l0| , and the numer- 
ical extrapolation of TBA results found for A £ 1R + , have led to the conclu- 
sion that the ground state energy of the SLYM has a square root singularity 
at i?(-A) 5 / 12 » 1.1325. In fact, this can also be seen directly from the initial 
TBA system (1.5). For this it is convenient to adsorb the bulk term into 



co(r), and work instead with the ground-state scaling function F$(r): 

2-7T r 2 1 

E (X,R) = -F (r) , F (r) = - -c (r) . (2.1) 



The function Fq thus defined is expected to be a regular function of r 12 / 5 fi"T[| , 
well-suited to analytic continuation. Negative values of A put r on the ray 
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r = pe 57 ™/ 12 , p € H + . A suitably-damped numerical iteration of the TBA 
equation (|1.5[) is convergent along this line, and the resulting ground state 
scaling function turns out to be real out to p = po ~ 2.99315, there being 
clear evidence for a square-root singularity at this point. This is shown in the 
lower set of points in figure |l], and matches in all aspects the TCSA results 
found earlier by Yurov and Zamolodchikov [flQf| . We conclude that the TBA 
is able to provide reliable (and, indeed, highly accurate) information even 
away from real values of r. 
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Figure 1: Two solutions to the basic TBA equation on the 
negative-A line (points) versus TCSA data (continuous lines) 

Next, the idea is to circle around the singularity, in the hope of picking 
up the next branch of the function - the first excited state. Analytic contin- 
uation is straightforwardly implemented when solving the TBA iteratively, 
by varying r step by step and at each new value of r taking the initial it- 
erate to be the final iterate of the previous step. The steps in r 
must not be too large, or the solution being tracked may be lost. Also, one 
must guard against numerical instabilities which are purely artifacts of the 
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iteration scheme. We adopted the simplest possible option, iterating as 



£ (n+i) = a J r cogh Q _ ^ L (n) + ( 1 _ ) E (n) (e) 5 (2.2) 

and finding empirically that values of a between 0.5 and 0.05 (depending on 
the value of r) gave optimal results. Stability rather than efficiency turns out 
to be the key issue, and it is worth seeking a more sophisticated approach. 
Nevertheless, fl2.2| ) was adequate for the current work. 

Of more profound import are the changes needed should a singularity of 
L{8) cross the real axis. In such a situation, further analytic continuation 
results in a function e which solves a modified TBA equation, with the 
contour of integration diverted away from from its original track long enough 
to avoid the singularity. To check for this, we must locate the complex zeroes 
and poles in z(9) = 1 + e~ £ ( 9 \ and be on our guard whenever any of them 
venture too close to the real axis. This is easily done, at least numerically - 
once e{8) is known along the real axis (or along some more general contour), 
equation ( |l.5|) provides an integral representation which can be used to 
reconstruct the function everywhere. The only point to watch is that the 
singularities of (j> at =tz7r/3 and beyond necessitate the introduction of extra 
terms when e(6) is continued beyond the strip —itt/3 < Im.9 < nr/3. 

Near the real axis at large values of |Re(0)|, the rcosh# term comes to 



dominate (1.5), and so two series of zeroes in z{8) are always seen, approach- 
ing the two half lines ±{6> G (D: Re# > 0, Im^ = [it/ 2 — Arg(r)]}. However 
if r lies on the initial segment of the negative- A line, Arg(r) = 57r/12 and 
< \r\ < po, we found evidence that, for the solution just discussed, a 
stronger result holds: up to our numerical accuracy, all of the zeroes in a 
strip along the real axis have imaginary part exactly equal to ±7r/12. As 
|r| = p — ► 0, the zeroes on the upper half line slide off towards +oo, while 
those on the lower head for — oo. The pattern becomes that of a pair of kink 
systems, one starting near 8 = — log(l/r), the other near 8 = +log(l/r). 
Conversely, as p — ► po, the two sets approach each other, although they 
always remain in their respective left and right halves of the complex plane. 

If r is continued in a clockwise sense about the critical point and back to 
the negative-A line, it turns out that none of these zeroes cross the real axis. 
The original TBA system ( |l.5| ) continues to hold, but its solution e(8) has 
nevertheless undergone a nontrivial monodromy. The values of the scaling 
function which result form the upper set of points in figure |], and match 
perfectly with the TCSA data for the first excited state, at least in the range 
of p for which our rather crude iteration scheme is stable. Hence the TBA 
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equation ( |1.5[ ) is at least doubly-degenerate along the negative-A line, with 
the second solution as physically relevant as the first. The monodromy has 
an interesting effect on the pattern of zeroes of z{0): while their imaginary 
parts apparently remain at ±7r/12, their real parts are no longer so simply 
arranged, at least for p smaller than about 2.8. The rightmost zero on the 
lower half line moves into the right half plane Re(#) > 0, thus lying to the 
right of the leftmost zero on the upper half line - the left and right kink sys- 
tems have become entwined, a feature that persists even as the two systems 
try to split apart in the p — > limit. On the other hand as p increases past 
2.8, the ordering is briefly restored, though there does not seem to be any 
great significance to this fact. Our numerics for the upper branch quickly 
become unstable in this region, but an extrapolation of the positions of the 
errant pair of zeroes is consistent with their moving continuously to the 
same positions ( ~ ±(0.24+z7r/12) ) as found for the first two zeroes on the 
lower branch, as p — > pq. This supports our supposition that the basic TBA 
system describes the first excited state on the whole segment < p < po, 
and not just on that part where our iterations converged. 

Having found the excited state for r on the negative-A line, we now 
continue r back to the real axis. As Arg(r) decreases from hit /VI, all of 
the zeroes of z{9) bar the first on the upper half line start to move up, 
towards the line Im.9 = ir/2. The first zero, on the other hand, is observed 
to move down, towards the real axis. The zeroes on the lower half line 
behave in a symmetrical fashion, as indeed they must given the 9 — > —9 
symmetry of the basic equations. This is precisely the situation mentioned 
above, and we should be ready to modify the TBA equations when the two 
singularities actually hit the axis. Unfortunately the iterative solution of 
the equation becomes unstable when singularities in L(9) get too close to 
the integration contour. It is possible to get a little further by distorting 
the contour along which the equations are being solved; in any event, a 
Pade extrapolation of the positions of the two singularities under suspicion 
clearly showed them crossing the real cLXIS ELS Arg(r) decreased. Assume 
that r is such that these two singularities in L(9) have crossed the real 
axis, and now lie at — 9q, 9q. Here and in analogous situations later on, 
we adopt the convention that of the pair {—9q,9q}, it is 9q which has the 
positive imaginary part after the axis is crossed. Coming in from the left, the 
integration contour for the convolution in ( |1.5[ ) must now first loop down and 
around the singularity at — 9q, and then back up and over the singularity 
at +6q, before proceeding to +oo along the real axis. An integration by 
parts turns these logarithmic singularities into simple poles; evaluating the 
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Figure 2: Pade extrapolation of 8q for the first excited state, 
from r = 1.5i (lower right) to r = 4.5 (upper left), and back. 



residues then allows the equation to be recast with the contour running 
along the real axis again: 



(2.3) 



e(9) = r cosh 9 + log lf a °°\ ~ <f>*L{6) 



Similarly, the expression for i*b(r) is modified, with cq of equation (|L 
being replaced by 



c(r) 



Ylr 



-i sinhf 



7T 



vr^ J- 



dOr cosh 6L(0) 



(2.4) 



Although these equations appear to contain an unknown parameter, namely 
#0) this is not the case: self-consistency demands that #o coincide with 
the position of the singularity in L(6) that necessitated its introduction in 
the first place. In this case, this translates as e(#o) = m (the branch of 
the logarithm to choose here can be fixed by continuity, starting from the 
situation before the singularity crosses the axis). Substituting 6 = 9q into 
(|2.3|) then gives 



= r cosh O - log S(29 ) - ^*L(9 ) . 



(2.5) 
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Flipping between equations ( |2.3[ ) and ( |2.5[ ), iterative schemes can be set 
up which are convergent in most regions of interest, given a reasonably 
accurate estimate of 6q for the initial iterate. Starting with the extrapolated 
singularities of the 'excited' solution to ( |1.5[) , a solution to (|2,4j ) can be picked 
up and followed all the way back to the real r axis. In this way, we tracked 
a solution along the line r = t + (1.5— t/3)t, with t varying from to 4.5. 



Equation (L5) converged for t less than about 0.45, while fl2,3j ) took over 
as soon as t became larger than 1. The problems for 0.45 < t < 1 seem to 
be artifacts of our iteration schemes, and in particular the good agreement 
the Pade-extrapolated singularity positions, both forwards and backwards, 
leave little doubt that the gap was crossed correctly. Figure § matches these 
extrapolations (shown by the continuous lines) with a selection of points 
from the 'raw' data. Points in the lower half plane (before the singularity has 
crossed the axis) derive from the unmodified TBA system ( |1.5| ), and those 
in the upper half plane from the modified system ( [2. 3D - (For the backwards 
fit, where a greater range of t was available along which to collect data from 
which to extrapolate, the matching is so good that the discrepancies with 
the target points are rather hard to spot on the figure.) It is also possible 
to monitor the behaviour of c(r), finding good agreement with TCSA data 
in both regimes. 

As r approaches the real axis, we see that 9q approaches z(7r/6 + <5(r)), 
with 5(r) a small positive correction which tends to zero for large real r. 
This information is enough to pin down the large-r asymptotics of ( |2.5| ). 
As r grows, the convolution term becomes relatively small, and the first 
two terms must cancel between themselves. The only way this can happen 
as r grows is for 28q to approach a singularity of S; with 9q ~ iir/6, the 
singularity at ot/3 is the relevant one. Substituting Bq(t) = i(tt/6 + 5(r)) 
and solving gives: 

o (r) ~ t (vr/6 + V3e-^ r / 2 ) . (2.6) 



This immediately tells us the leading asymptotic of c(r): substituting (|2.6|) 
into first term in ( [2.4D gives c(r) ~ — 6r(l + 3e~ v ^ r / 2 )/ir. The next term is 



also easy to find, dropping the third factor in (2.3) and using the zeroth- 
order value of #0; namely iir/6, in order to find the leading behaviour of e{9) 
for substitution into (^J|). Gathering everything together we find 



(2.7) 



c ( r ) —11 + 3e-^ r/2 / d0 coshf? S{6 + i7r/2)e~ rcosh6 
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where the bootstrap relation S(9 + m/6)/ S(9 — m/6) = S(9 + m/2) , relevant 
because 9q has been given its asymptotic value, was used to reduce the 
two S-matrices in (|2.3|) t o one. This matches exactly with the asymptotics 
predicted in refs |10|, |12| for the spin-zero one-particle state, and therefore 
lends strong support to our proposal. Further evidence will come from the 
numerical comparisons with TCSA data to be reported shortly, but first 



we would like to mention a natural generalisation of equations (2.3) - (|2.5| ) 
which seems to capture all the remaining one-particle states. 
The equations to consider read, for r 6 H, as follows: 

s(9) = r cosh 9 + log S Jf~l°\ -<f>*L(e) ; 

6r — 3 r°° 

c(r) = i— (sinh0 o -sinh0 o ) + -o / d9 r cosh 9 L( 9) , (2.8) 

7T 7T Z J-oo 



where 9q, 9q are the complex-conjugate locations of a pair of singularities 
in L(9). We will discuss the one-particle states with positive spin, and so 
we take Re(#o) > 0; the negative-spin states work similarly. The earlier 
equations at real values of r are recovered if 9q is forced to be purely imagi- 
nary; conversely if r were to become complex in (^^), then 9q and 9q would 
generally cease to be complex conjugates, and would have to be tracked indi- 
vidually. Equation (|2.5| ) must also be modified, both because of the changed 
form of ( |2.8| ) and also to allow for a more general singularity at 9q, namely 
s(9q) = (2n+l)iri (and n > for Re(#o) > 0). The equation becomes: 

2niri = r cosh 9 - log S(2ilm(9 )) - <p*L(9 ) . (2.9) 

Part of the motivation for these equations came from the form of the 
two-particle equations to be introduced in the next section; but to see im- 
mediately that they have a chance of being correct, consider the large-r 
asymptotics. The convolution term is sub- leading, and again it turns out 
that the balance between the first two terms is achieved via Im#o = n/6+5(r) 
with 6(r) vanishing as r — > oo. Now take the imaginary part of ( |2.9|) (previ- 
ously this vanished automatically for real r), and consider its behaviour as 
r becomes large: 

2nvr = rsinh(Re(6» ))sin(7r/6+5(r)) - ImlogS(2?(7r/6 + <5(r))) 

~ ^smh(Re(0 o ))-|(l-sign(5(r)))7r . (2.10) 
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The term ^(1— sign(5(r)))7r is included to allow for 6(r) being negative, in 
which case S(2z(7r/6 + S(r))) is negative and the logarithm picks up an 
imaginary part. Feeding this into (^j^) gives 

6r 

' ' ' (2.11) 



or i 

c(r) Jl + {2-Ks/r) 2 

TT v 



with s = 2n + i(l— sign(<5(r))). This is exactly as expected for a one-particle 
state with spin s. 

One final observation: for all but the spin-zero one-particle state, equa- 



tion (2.8) is not symmetrical under 9 — > —0. Hence it can never be obtained 
by analytic continuation of the ground-state equation. But this is just as one 
would expect: the entire Hilbert space of the SLYM splits up into sectors of 
different spin, and analytic continuation can only ever move levels around 
within a given sector. 
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Figure 3: Proposed one-particle scaling functions (points) com- 
pared with TCSA data (continuous lines) 

Figure || compares the numerical solutions of equations (^) , ( pT9| ) with 
TCSA data, for the first four one-particle levels. We used the TCSA program 
of ref. fl3[| , truncating the quasiprimary fields at level 5 and including as 
many derivative states as the program allowed. The discrepancies are too 
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small to see on the figure, being of order 10 at small values of r, growing to 
10~ 3 — 10~ 4 (depending on the level of the state) at r=10. The former errors 
can be traced to our numerics and should not be too hard to reduce, whilst 
the latter can be ascribed to truncation effects in the TCSA - in particular, 
they grow as the level of the state under consideration gets higher. 

The only problem that the reader might notice is that the lowest set of 
points stops short, at r 2.53. By this stage S(r) has become so large that 
#o has almost reached in/ '3. We should expect trouble at such a point as 
a new kind of singularity enters into the equations, caused by singularities 
in S(9—8q) and S(0 — Oq). We have tried to take this into account, but 
the resulting equation appears to be much less stable and has resisted our 
attempts at an iterative solution. We will return to this problem in section 4; 
in any event it only appears to trouble the spin-zero state. 

3 Two-particle states and beyond 

The singularity at ro = poe 5m ^ 12 is only the first of a whole sequence of 
singularities seen by the 'zero-particle' TBA ( [1.5|J1.9| ). They are approxi- 
mately evenly-spaced along the direction of the imaginary axis: the next is 
at ri Ri 0.5311+9. 1346?, and the next at r 2 « 0. 42+15. 44?. Figure f| shows 
the positions of ro and r±, as emerged from the numerical solution of the 
basic TBA equations on a suitably-fine grid. 



10 8 6 4 2 




Figure 4: A contour plot of lm(F(r)) as obtained from the basic 
TBA equations (1.5), (1.9), with r lying in the box < Re(r) < 1.85, 
< Im(r) < 10.1, showing the branch points at ro and r±, and also 
the segment of the negative- A line, running from to ro, along which 
Im(.F(r)) = 0. The box was scanned from right to left, so it is the 
lower branch of this segment that is visible. 
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On the analytical side, we note that the first iterative correction to the 
Ising-like behaviour of (|1.5| ), namely 4>* log(l+e _rcosh6 '), tends uniformly to 
zero as Im(r) — > ±00 with Re(r) > held fixedQ Hence we expect that 
the locations of the r n will eventually approach (2n+l)-7rz as n — > 00. These 
appear to be the only branch points on that part of the Riemann surface of 
F(r) explored by the zero-particle TBA. However the full surface must have 
much more structure: the behaviour of the action (|1.4| ) under A — > A means 
that each singularity r n beyond the symmetrically-placed ro must have an 
image on some other sheet, located at 

r n ee e 5 "/ 6 r n . (3.1) 

These images, invisible to the zero-particle TBA, should instead be seen by 
the more general TBA systems that we are trying to construct. 

This idea is confirmed by the values of F{r) which result from the one- 
particle TBA introduced in the last section: a grid plot of its real and 
imaginary parts exhibits a square root singularity at f\ rj 4. 1074+8. 1763«. 
An example of such a plot is shown in figure ||. (Incidentally, that this works 
is a further piece of support for the results of the last section.) Judicious 
use of the TCSA, allied with TBA data as a check on accuracy, is a great 
help in mapping out how the various sheets fit together, and seems to be 
qualitatively reliable at least out to |r|=20. Here, it tells us that the branch 
point at f 1 should connect with the first two-particle state. We can therefore 
play the same game as before, and study the behaviour of the solution e{9) 
of the one-particle TBA as r is continued round r±. This should yield a 
two-particle TBA equation. 

When r approaches the branch point, we found that the singularities 
at 9 = =t#o, already implicated in the one-particle equation, remained near 
±z7r/6, whilst a second pair headed for the real axis, just as happened when 
the one-particle TBA was being formed. Strictly speaking we should now 
repeat the rest of the one-particle work, this time keeping track of two 
independent singularity positions, #0 and as r returns to the real axis. 
Numerically this is delicate: with both 9q and 9\ to locate, an efficient 
iteration scheme is hard to find, and so we will leave this question to one 
side for the time being. Besides, in some respects the key piece of information 
has already been obtained: the two-particle TBA equation should involve 
four singularity terms, tied to singularities in L{9) at ±9q and ±9\. Once r 
reaches the real axis, the situation changes favourably: we expect the four 

1 when checking this, it is helpful to note that log(l+ae ,7 )d7 = for \a\ < 1 . 
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Figure 5: Im(F(r)) in the box 3.8 < Re(r) < 4.4, 7.7 < Im(r) < 8.5 
as obtained from the one-particle TBA equations (2.3), (2.4), showing 
the branch point at f\. 



singularities to be invariant not only under Oi — ► —6i, but also under 9i ^ 9{. 
This reduces to one the number of independent singularity positions, making 
a numerical solution no harder than the cases examined in the last section. 
The TBA equation to solve reads: 

m = rcosh9+log |^ + log §|±M-,, iW; 

c(r) = i (sinh# -sinh# ) + — / dQ r cosh 6L(0) , (3.2) 

7T 1T Z J-oo 



with #o satisfying 

e(0 o ) = (2n + 1)ot (3.3) 

For the lowest two-particle state, found on continuing round f%, we found 
e(9q) = nr, so that n is equal to zero. However, the results of the last section 
make it very natural to allow for the more general possibility, in the hope 
of catching the other two-particle states. As with the one-particle TBA 
equation ( |2.8j ), #o and #o cease to be complex conjugate if r strays from the 
real axis, eventually metamophosing into #o and —B\ as f\ is approached. 

The analysis of the infrared limit proceeds in much the same way as for 
the s ^ one-particle states. Substituting 6 = 8q in the first of ( |3.2j ) and 
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taking the imaginary part of the large-r asymptotic gives us 



27r(2n + |(l-sign(g(r)))) ~ r sinh (Re(fl )) + 2Imlog v " u; (3.4) 

where, as before, Im0o = vr/6+5(r) and 6(r) — > as r — > oo. In the functions 
5, can be replaced by its limiting value, so these terms become 

ImlogS(2Re(0 o ) +»tt/3) = 5lmlogS(2Re(0 o )) (3.5) 

(using the 3 property of 5(0) ) and 

Im log S(2Re(0 o )) = -« log S(2Re(0 o )) (3.6) 

(Recall that S(0) is a pure phase for 6 real, so the right-hand side of ( |3.6| ) 
is indeed real.) Combining all of these terms together, ( |3,4| ) becomes 

2^(2n + ±(l-sign(<f (r)))) ~ r sinh (Re(0 o )) - t log S(2Re(0 o )) • (3.7) 

This is just the Bethe-Ansatz quantization condition for a two-particle state 
with rapidities (— Re(#o)> Re(#o)) an d Bethe quantum numbers 

(-2„ + I si|! n(i), 2„-Isign W ) = (-i, i) , (-§. | ) , (-§, jj ) . . . (3.8) 

with n = 0, 1, 1, 2, 2, . . . and sign (5) = —1, 1, —1, .... (cf, for example, equa- 
tion (4.5) of ref. pQ| ). These values of sign(£) should be imposed when 
handling the equation numerically, so that when seeking to follow a particu- 
lar solution, the idea is to specify not only the value of e(0o)> equation (|3.3|) , 
but also the sign of S(r). The only point where we were not able to impose 
a particular sign on 5 and obtain reasonable results was when we attempted 
to set sign(<5) = +1 when n in ( |3.3| ) is equal to zero. From a physical point 
of view, this is as it should be, given the exclusion principle - otherwise, the 
state (— |, |) would have appeared twice. 

Substituting the value of Re($o) which follows from (|3,g|), together with 
Im(#o) = 7r/6, into ( |3.2j ) automatically gives the expected algebraic asymp- 
totic for F(r). The next section will examine the ultraviolet limits analyt- 
ically; in the meantime figure |6| illustrates the numerical agreement with 
TCSA data for the first three states. The discrepancies were much the same 
as for the one-particle sector, described earlier. 
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Figure 6: Proposed zero-momentum two-particle scaling func- 
tions (points) compared with TCSA data (continuous lines) 



By this point the structure is clear enough that we can conjecture a 
general equation for an n-particle state: 

e{0) = r cosh + ^ log S Jf ~ y - foL(fl) ; 

i b(U — t>i) 

6r — 3 f°° 

c (r) = »V— (sinh0i-sinh0i) + -5 / d0rcosh0L(0) , (3.9) 

We expect that particular states are selected by imposing the values of e(0i), 
and also the signs of the various Si. We only analysed in detail the spin zero 
three-particle states, where steps similar to those already described give the 
following Bethe Ansatz quantum numbers: 

(-2n + \ + isign(<5), 0, 2n - \ - |sign(5)) = (-1,0, 1), (-2, 0, 2), . . . (3.10) 

with n = 1, 1, 2, 2, 3... and sign(<5) = +1, —1, +1, .... Again, the exclusion 
principle leads us to suppose that the n=0 case should be omitted. 

This large-r asymptotic, together with the special cases already studied 
in some detail, lend support to ( |3.9| ). However numerical work will also be 
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needed. In particular, the equation will have to be modified at values of 
r below some critical r c whenever the phenomenon observed above for the 
spin-zero one-particle state occurs, and the imaginary part of a singularity 
position Oi ventures too far from ir/6. We suspect that this will happen 
whenever one of the Bethe Ansatz quantum numbers is equal to zero, and 
thus will afflict the three-particle states just discussed. In fact there are signs 
that this is just what is needed to get the ultraviolet asymptotics correct in 
these cases, a point that we shall return to later. 



4 Ultraviolet behaviour 

Analytic work has thus far been restricted to the infrared regime. In this sec- 
tion the opposite limit is considered, and the ultraviolet scaling dimensions 
of various states extracted. We start with the two-particle states. 

In the limit r — ► 0, the two-particle TBA splits into a pair of kink 
systems, just as for the ground state. For the right kink system, the terms 
from the singularities at —9q and —9q drop out, and r cosh# can be replaced 
by \re b . The real part of 9q tends to infinity like log(l/r); we will also need 
a little information about the behaviour of the imaginary part. Recall that 
for large r, Im(#o) ~ ir/6+5(r) with S(r) — > as r — > oo. In the opposite 
direction, we find that, although 5(r) grows, it tends to a finite and still- 
small limit as r — > 0. The precise value drops out of the final equations; 
the only important fact is that for all of the two-particle states (and also for 
all of the s^O one-particle states) its absolute value is less than tt/6. Quite 
why this is important should become clear shortly. The scale-invariant kink 
TBA reads: 

e(e) = r -e s + lo g fl^- ( p*L(e). (4.1) 

(Note, to pass from this 'chiral' TBA equation to the general one-particle 
equation introduced in section 2, we just need to replace with cosh#.) 
Next, eliminate r by replacing 9 with 9 — logr and e(9) by e(9 — logr), to 
find 

e(9) = \e e + log - ^L{9) , (4.2) 

where rj = #o+logr. In terms of these quantities, the limiting value of c(r) 
is 

c=4/ d9e e L{9) + -i{e r > -e*!) . (4.3) 
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Now we proceed in the usual manner. First take a derivative with respect 
to 0: 



89 w 2 89 * S(9 -rj) 
Now substitute this into the formula for c: 



50 



Then (remembering that 0(0) 



-zj| log 5(0)) 



(4.5) 



6 



7T 



(fx log(l + e x ) + -e min log(l + e" 



2cf)*L(r]) - 2(p*L(7]) - e* 1 + e v 



(4.6) 



The first piece can be expressed in terms of the Rogers' dilogarithm function 
C(z) = -lf*dt(^ + ^) as 



7T- 



£ 



1 + e £ 



1 



whilst the second can be evaluated using 



1 + e Emc 

for 77 and ?y: 



(4.7) 



2niti 
-2niri 



log(|5(z2Im(r/))|) - 7ni(l-sign(5)) - 0*L(t?) 



e* 7 + log(|S(*2Im(r/))|) + vu±(l-sign(<5)) - </>*L(r?) , (4.8) 



and is equal to 12(4n+l— sign (J)) . 

Finally, e m in is reached as 9 — > — oo, e ma x as — > oo. Since S(±oo) = 1, 
the limiting form of (4.1) gives £ m i n = log(l+v5)/2) and e max = oo, 
a result which was also checked against our numerical solutions. Since 
£(2/(3+ v / 5)) = vr 2 /15, this gives us 



12(4ra + 1 - sign(5)) 



(4.9) 



The calculation for the s ^ one-particle states is essentially identical. The 
only difference is that for one of the two kink systems is in its 'ground state', 
in that there are none of the extra terms involving 9$. Hence the size of the 
additional contribution is halved, and 



6(4n + 1 - sign(<5)) 



(4.10) 



17 



Recalling from the infrared result ( |2.11 ) that s = 2n+^(l— sign(<5)), this is 



-1/12 times the scaling dimension on the cylinder of a spin s descendant of 
the primary field ip, as expected. 

For s = in the one-particle sector, the situation changes and our discus- 
sion must be much more tentative, as our numerical work lost stability below 
r = r c rs 2.53. Nevertheless we can conjecture a plausible modification to 
the equation which seems to predict the correct ultraviolet asymptotic, and 
the remainder of this section is devoted to this question. 

First, we must decide just how we expect the equation to change. In par- 
ticular is important to know if the singularities at ±9q stay on the imaginary 
axis, or if they somehow each split in two and acquire a real part. Unlikely 
as it appears at first sight, it seems that it is the latter possibility which ac- 
tually occurs. One can gain information about this by again returning r into 
the complex plane, and following the position of 9q as r returns to points on 
the real axis below the critical value r c . Pade extrapolation clearly showed 
a finite limiting value of Re(#o) in each case we tried, the value increasing as 
the limiting r decreased. For example, extrapolating down from the upper 
half plane, we obtained Re(6>o) 0.77 at r=1.5, and Re(#o) ~ 1-17 at r=l. 
The extrapolation of Im(#o) was not so reliable, but the results were consis- 
tent with the limiting value being ir/3 every time. In any event, this seems 
to be the only way to make sense of the non-zero limit for the real parts, 
and so we will proceed on the assumption that indeed lm(9o(r)) = tt/3 for 
all r < r c . 

When Im(#o) finally reaches vr/3, the singularities in L{9) caused by the 
log S{6 ± #o) term hit the real axis, forcing a deformation of the contour of 
integration. Integrating by parts in the convolution, and then taking the 
principal part to restore the contour to its original track, equations pi 



(2.4) should be modified thus: 



e(0) = ^osh^log^^ + llog ^^^^^ -4»W) ; 



12/* Qt 
c(r) = i sinh #o — % — sinh(#o — in/3) 

TT TT 

+ — d9 r cosh 6 L{9) , (4.11) 

7T J — oo 



The factors of 1/2 appear because the singularities remain on the real axis. 
This equation looks rather complicated, but simplifies once we put 9q = 
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/3, with (3 real. Using the bootstrap equation we find 



,. S(9-/3-m/3) !, S(6+f3-iir/3) , r ,„, 



c(r) = ^/3cosh/3 + — / dQ r cosh9L(9) . (4.12) 

7T 7H J-oo 



To understand what has happened, it is worth thinking about how the 'sin- 
gular values' of e{9) have moved around. This can be discussed using the 
functions Y(9) = e £ ^ which solve, for the Lee- Yang model, the following 



functional equation or Y-system [11] 



Y{9-i^)Y{9+i^) = l + Y{9) . (4.13) 

The special points that we are interested in, 9q being an example, are those 
for which 1+Y = 0. Given such a point 9^°\ the Y-system can be used to 
find a sequence of other points 9^ = 9^+km/3 where Y also takes special 
values. Setting Y^ = Y(9^), with y( fc + 5 ) = y( fc ) by the periodicity of 
the Y-system, the sequence cycles round as —1,A,B,—1,0, with A+B = 
and Y^ ) either the first or the fourth term. For #0j the former option is 
realised and so as 9q moves around, it carries with it the points 9o—nr/3 
and 9q— 2^7^/3, at which Y takes the values and —1 respectively. While 
r > r c , these points lie on the imaginary axis, together with a symmetrically- 
placed triplet of points associated with — 9q. When r reaches r c , the points 
9q— 2z7r/3 and —9q collide at —itt/3. Since for both of these points 1+Y = 0, 
this opens up another possibility to distort the singularity positions while 
maintaining the symmetry under 9^9 that they must respect while r 
remains real. This seems to be realised here: as r decreases below r c , the 
points #o — 2z7r/3 and —9q move away parallel to the real axis, with a similar 
story at +iir/3 ensuring that the symmetry is preserved. By continuing 
r away through complex values, we were able to observe this end result 
without having to tangle with the particularly singular behaviour actually 
at r=r c . 

The last task is to find an equation for 9q, or equivalently for j3. The 
value of e(9o) remains equal to in, but one must be careful when sub- 
stituting 9=9q into the first of ( 4.12j ). This is because the singularities 



just discussed cause the right hand side of this equation to develop a cou- 
ple of singularities precisely when 9=9q. Of course, the overall result re- 
mains finite, and the simplest approach seems to be to consider the limit 
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of e((3+iit /3—ie) — e(f3—nr/3+ie) as e -> + . From one point of view this 
is equal to 2^7^; equating this with the result from the right hand sides we 
found 

m = iy/Er SmhP + 2 bg 5(2/? + Z2VT/3) S(20 - ,2ir/3) 

- (0*L(^ + «r/3) - 0*L(/3 - ct/3)) . (4.14) 

(We should mention here that we have been somewhat cavalier throughout 
in our treatment of the branch choices for the logarithms. This issue deserves 
a careful study, especially in regard to the way the branches behave under 
analytic continuation between energy levels.) 

This concludes the modifications to the one-particle equation. If the 
comment made at the end of section 3 is correct, then a similar manoeuvre 
will be needed whenever a zero-momentum particle is present in a Bethe 
Ansatz state. The natural generalisation of the equations just obtained for 
such situations, correcting equation ( |3.9[ ) for r less than some r c , is 

W + 2 g S{9-(5 + zvr/3) S(6 + p + ztt/3) 

c(r) = V^cosh/? 

7T 

12r — 3 f°° 

+ %Y^ (sinhfli-sinhfli) + — / d6rcosh6L(0) , (4.15) 

7T 7T J— oo 

where as above the singularity position 9q ceased to be purely imaginary at 
r=r c , and was replaced in the equations by the real variable (3 = 9q—17t/3. 

Now we must analyse the r —* limit. The first thing that we notice is 
that the modification has achieved the remarkable trick of splitting the 'zero- 
momentum' singularity position 9q into constituent singularities at z7r/3±/3, 
which can now join the respective left and right kink systems. For the one- 
particle case, assuming that (5 — ► oo as r — > 0, and replacing 9 with 6 — logr 
and e(9) with e(9 — logr) , the right hand system becomes 

Q Q rod 

c = --VSe^ + A, d9e d L(9), (4.16) 

7T 7T Z J-oo 



20 



where r\ = (5 + log r satisfies the 'quantization condition' 



V3 



itt = % e ' 



!>*L(t7 + z7r/3) — <f>*L(r] — z7r/3)) 



(4.17) 



Now the calculation runs as before, modulo one subtlety to be mentioned 
shortly. One finds 



3\/3 _ 6 

6 ~t~ n 



TT 

6 

TT 



7T 



dx log(l + e x ) + Ismin log(l + e~ 



? [0*Zy(77+z7r/3) — (f)*L(rj— wr/3)] 



(4.18) 



Recognising the dilogarithm and using (4.17), 



6 



7T^ 



c 



1 



l+e £ 



1 



As before, e max = oo, and e £min solves e 2emin 
has two solutions: 

1 ± \/5 



-6. (4.19) 
an equation which 
(4.20) 



Previously we selected the positive solution, e{9) being a monotonically 
increasing real function in that case. However this time note that e^ 00 -* = 
e £max — anc j that as # decreases along the real axis e £ ^ falls all the way 
down to zero at 9 = (3. This suggests that the relevant solution by the time 
9 = — oo is reached is the negative one (note, such solutions were previously 
observed to be of relevance to the excited states in some other models by 
Martins, in ref. ||]). The correct prescription for C(x) for general x £ H + 
is (see for example ]l4| ) 



C(x) 

which in this case gives 



TT 



C(l/x) for x > 1 , 



TT 2 4-7T 2 

£(2/(3-^5)) = - - C(2/(3+VE)) = — 



(4.21) 



(4.22) 



This gives us c = —22/5, as expected for the spin-zero one-particle state, 
which in the ultraviolet is created by the identity operator /. (Note, the 
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theory being non-unitary, this is not the ground state!) We also analysed 
the spin zero 3-particle sector, and found 

22 

c = - — - 12(4n+l-sign(<5)) (4.23) 
o 

with n = 1, 1, 2, 2, . . . and sign(e)) = +1, —1, +1, . . . just as in the infrared. 
Note that here also n must start from 1: the two options at n=0 give either 
—22/5, double-counting the state /, or —22/5 — 24, the would-be scaling 
dimension on the cylinder of the null field ddl. This seems to suggest 
a relation between the exclusion principle on the Bethe Ansatz quantum 
numbers, and the null field structure of the conformal states. 



5 Conclusions 

This work is still in its early stages. We would like to have a more secure 
understanding of the zero-momentum one-particle state, and in particular to 
be able to follow its behaviour numerically all the way down to r=0. Work 
on this question is in progress. The situation for the remaining one-particle 
states, and for all of the two-particle states, is very satisfactory but beyond 
that our equations become more conjectural, albeit natural. Numerical and 
analytic work is needed, both to confirm their status and to unravel their 
structure. On the numerical side this poses the particular challenge of devel- 
oping efficient methods for the tracking of a number of singularities simul- 
taneously, but this should not be insurmountable - the three-particle state 
would be a good starting point. Quite apart from the analytic insights we 
can hope for, the method promises to be very competitive numerically with 
the TCSA, particularly for higher levels and larger values of r. Note also 
that all of the results we have obtained for the SLYM are directly relevant, 
after a multiplication by two, to the thermally-perturbed three-state Potts 
model. This follows from the simple relationship between their respective 
ground-state TBA systems Q. 

Turning to more general issues, many new features of the TBA equations 
seem to emerge when the whole complex r plane is considered, and there 
remains much to explore. The map A — > A, simple in its effects from the 
point of view of the perturbative action ( |1.4| ), is far from trivial when acting 



through equation (3.1), as r —* f, on the space of multiparticle TBA equa- 
tions. More locally in this space, the study of the zero, one and two particle 
TBA systems has shown how one equation can melt into another with the 
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passage of singularities across the real axis, and it is important to extend 
this treatment further. An immediate hope would be to thereby justify the 
more general equations ( |3.9| ) conjectured above. More ambitiously, given 
the good understanding of the ultraviolet and infrared limits of the model 
provided by conformal field theory and the Bethe-Ansatz quantization con- 
ditions respectively, one might try to say something about the structure of 
the full Riemann surfaces for F(r) in the various sectors, and to understand 
the way in which the domains of the various TBA equations are patched 
together on these surfaces. 

We would like to stress how general the method advocated in this paper 
should be. As a first step, there seems to be no serious obstacle to its 
application to the known purely elastic scattering theories. In the absence 
of any numerical work to report at this stage, we will restrict ourselves to 
a couple of comments. First, we note the crucial role that the 4> s property 
of S(0) played in much of the analysis presented above. One can anticipate 
that more general cases will exhibit a similar interplay between the algebraic 
properties of the S'-matrix and the asymptotics of the multiparticle TBA 
equations. Second, many scattering theories, purely elastic and otherwise, 
exhibit additional symmetries which on the one hand divide the Hilbert 
space into further subsectors, and on the other permit the construction of 
alternative 'seed' TBA Qj. These will provide additional starting-points for 
the continuation process, thus helping to fill out the extra sectors. In the 
Ising model example discussed in section 1, this was exactly how the states 
related to the spin field arose. 

Our main conclusion, however, is rather more general than this. It is 
that any ground-state TBA equation encodes within itself equations for 
many excited states, and that analytic continuation provides the means by 
which these equations can be extracted. It will be interesting to see just 
how far this programme can be carried through. 
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Note - Ref. [15], by Bazhanov et al, appeared as we were finishing this 
paper. These authors independently obtained TBA-like equations for the 
spin- zero excited states of the SLYM, though by a very different route from 
the one that we have been describing here. 
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